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Abstract We study the boundary behaviour of solutions u of — Ajvit + 1u|'~^m = in a bounded 
smooth domain f2 C subject to the boundary condition u — Q except at one point, in the range 
q > N — 1. We prove that if g > 2N — 1 such a u is identically zero, while, if — 1 < g < 27V — 1, 
u inherits a boundary behaviour which either corresponds to a weak singularity, or to a strong 
singularity. Such singularities are effectively constructed. 



1 Introduction 

Let ri be a domain is {N > 2) with a compact boundary d^. Let 5 be a continous real 
valued function and a G dil. This paper deals with the study of solutions u € C^{tt \ {a}) 
of the problem 

-div(^\Duf^^Du^+g{u) = inQ 

u = on dft \ {a}, 

and we shall be more specifically interested in the case when g has a power growth at 
infinity. When N — 2, this problem falls into the scope of the boundary singularity problem 
for semilinear elliptic equations. The study of the A'^-dimensional problem 

J -Am + g{u) = in fi , , 

\ u = ondn\{a}, ^ ' 

has been initiated by Gmira and Veron in [7]. Among the subjects under consideration 
were the question of removability of isolated boundary singularities and, in the case such 
singularities do exist, their precise description. This seminal article was at the origin of a long 
series of further works by Dynkin, Kuznetsov, Le Gall, Marcus and Veron in the framework 
of the trace theory and, later on, the fine trace theory in the case where g{r) = r \r\'^~^, 

*To appear in Journal of Functional Analysis 
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g > 1. One of the main reasons for such a large impact consists of the observation of the 
existence of a critical exponent q = q* = {N + 1)/{N — 1). If g > g* any solution of 

-Am + ^ u = Q in il n -j^ 

u^O on 917 \ {a}, ^ ^ 

is identically zero, while if 1 < g < g* it appears that there exist two possible behaviours of 
singular solutions near a, the solutions with weak singularities and the ones with the strong 
singular behaviour. Later on, these two types of singular solutions played a fundamental 
role in the description of the rough trace of positive solutions of p. 3 |) . 

Although the techniques needed are considerably more refined, it appeared that the 
description of solutions of (jl.l p inherits the same structure as for (|1.2 p . The first step is 
to understand the model case problem 



(1.4) 



(1.5) 



div (\Du\^ ^ Du^ + \u\'^~^ u ^ infl 

u = Q on dfl \ {a}, 

To this equation, we associate the homogeneous equation 

- div (\Duf^^ DuJ =0 inn 

u = on dfl \ {a}. 

It is proved in [5] that for any fc > there exists a unique solution u = Uk oi (jl.5 p satisfying 

Uk{x) = k ^^^^ ^ (1 + o(l)) as 2: — > a, (x — a)/ \x — a\ ^ a, (1.6) 
\x — a\ 

where p{x) — dist {x,dil). When fc = 1, this solution plays the role of the Poisson kernel, 
although neither any weak formulation nor any reasonable trace theory seems to exists, and 
we shall denote it by V^. The behaviour (|1.6 p (up to a multiplicative constant) corresponds 
to weak singularity behaviour for (jl.l p . whenever such singularities exist. The first result 
we prove is the following: 

Theorem Let N — 1 < q < 2N — 1 := gc • Then for any k > there exists a unique solution 
u — Uk^a of problem {1.^ [ ) satisfying 

Uk.aix) = k ^ ^ {1 + 0(1)) as X ^ a, (x — a)/ |x — a| — > a. (1-7) 
|a; — a\ 

Furthermore Uoo,a = linife_,oo exists and is a solution of {H [ ) which satisfies 

hm \x-a\''''^'^+^-''K^^a{x)=uo{<j), (1.8) 



X — > a 
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and to is the unique positive solution of the following quasilinear equation on the upper 
hemisphere of the unit sphere S^~^, 



- A (^(3^y + 



2\^^-'^/\, 19-1 n Q^-i (1-9) 
' a;+|Lj| uj = on 

w = on dS!^-\ 



where Pq — N / + 1 — N) and A = {N — 1)0^. The proof of the existence of Uk^a, as well as 
its singular behaviour, is settled upon the conformal invariance of the A^-harmonic operator 
and the construction of subsolution of the same equation. Estimate (jl.8 \ is proved by 
scaling method. The role of the critical exponent q^ = 2N — 1 is enlighted by the following 
result. 

Theorem Let g he a continuous function such that 

(i) liminfr^oo > 

(1.10) 

(m) lim sup^^_^. g(r)/ Ir]''' < 0. 
Then any function u £ C^{Q, \ {a}) solution of il.l \) extends as a function u G C(f2). 

As in the semilinear case, the occurrence coincides with the case where the blow-up 
exponent —f3q which is natural for equation (|1.4 p coincides with the one of the function 
solution of (|1.5 p . Finally we provide the full classification of positive solutions of problem 

(m. 

Theorem Let N ~ 1 < q < qc and u is any nonnegative solution of \H [ ), then 
(i) Either u = Q, 

(a) Either there exists k > such that u = Uk.a- 
(Hi) Or u ^ Uoo,a- 

In the proof of (iii) the boundary Harnack inequalities that satisfies any positive solution 
of (|1.4 p (see [2]) play a fundamental role. 

Our paper is organized as follows 

1- Introduction 

2- Weak and strong boundary singularities 

3- The removability result 

4- The classification theorem 

2 Weak and strong boundary singularities 

The construction of positive solutions of 

- div (\Duf-^ Du^ + u = 0, (2.1) 
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is settled upon three facts: the existence of sohitions to the homogeneous equation 

-div(\Duf^^Du^^O, (2.2) 

the conformal invariance of (|2.2 p and an a priori estimate satisfied by any solution of (|2.1 p . 
Throughout this paper C denotes a positive constant which depends only on the structural 
assumptions corresponding to N, p, q and fJ. The value of the constant may change from 
one occurrence to another. 

Proposition 2.1 Let Vl CZ M.^ be a domain with a compact boundary and a € 917. Consider 
real numbers q > p ~ I > 0, A > and B > 0. If u e C(Jl \ {a}) n Wl;^{n) is a weak 
solution of 

i -div (\Duf^^ Du^ + A\u\''^^ u < B mfl 



u < on dfl \ {a}, 



(2.3) 



it satisfies 



^i^) < [ ^TT—Jp) +(^) Vxel7\{a}, (2.4) 



where A and /i depends on N , p and q. 
Proof. By assumption 

/ (\Duf~^Du.D( + A\u\'^~'^uC)dx<Bf(dx (2.5) 

for any ( € W^'P{^1) with compact support, C > 0. Let rj £ C^(M) be a nonnegative function 
such that < 77' < 1, 77" > 0, 77 = 77' = 77" on (-00,0], < r]{r) < r on (0,cxd). For 
e > we set rjeir) = rj{{r — e)+). Let C, G M^^'P(R^ \ {0}) with compact support. Inasmuch 
{ri'^(u)Y^^C, has compact support in Q, and 

D {{v'M)Y-'0 = {v'M)Y-'DC + {p- iMiu)r-'i:iu)CDu, 

it belongs to W^'P{n) and is an admissible test function for (|2.5 |) . Thus 

J {\Duf-^ Du.D {iv'Au)y~^C) + A |7i|«-^ u{ri{u))P-^C^ dx < B j {rj'^{u))P-\dx , 

and 

iDul"-^ Du.D {{r^'^{u)Y~\) > {r^'^{u))P-^ \Du\"-^ Du.DC = \Dv,\''^^ Dv,.DC, 

where we have set — r]^{u). Furthermore, r] can be chosen such that r'^ {rj'^{r)Y~'^ > rij{r), 
for example if we fix ri{r) — r^/26 on {0,6] and ri{r) = r — S/2 on [(5, 00) for some (5 > 0. 
We extend by outside \ {a} and denote by Ve the new function, then e W^/oc (M^ \ 
{a}) n C(R^ \ {a}) and 

J (\Di^f~^ Di^.D(^ + A\i,\'^'^ i,C)dx <B j Cdx. (2.6) 
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This means that -Dg is a weak subsohition in \ {a}. By [181 Lemma 1.3], we derive 

for some A > and /i > depending on A'', p and q. Letting successively e — > and S ^ 
we obtain (|2X| . □ 

When Q is smooth we have a sharper estimate 

Proposition 2.2 Let 17 C be a bounded domain with boundary and a £ dil. Let 
q > p - 1 > I and a > 0. If u € C(U \ {a}) n Wl^^{n) is a weak solution of k2.3 \) with 
B = 0, there exists C > depending on fl, p and q such that 

n{x)<- ^^^':\/(^^,^^) Vxen\{a}, (2.7) 



where p{x) = dist (x, d^l. 

Proof. By translation we can assume that a — Q. For e > let be the solution of 



-div 



= w+ on dVL'^ . 



By [181 Lemma 1.3] as in the proof of Proposition 12.11 and the maximum principle, there 
holds 

Consequently t < t' Ve > v^i. Letting e — > and using the previous inequalities and 
the classical regularity results for solutions of quasilinear equations [12] we conclude that 
converges, as e ^ 0, to some v which is a nonnegative solution of 

f -div {\Dv\^'^'^ Dv^ + Av'i = 0, in n 

\ v = Q ondVL\{Q}. 

and dominate u. Further, if £ > the function defined by v^{y) — iP^^i^^~P'>v{iy) is a 
solution of (|2.9 p with ft replaced by il^ = Let x e fl \ {0} and i = \x\. Since 

)^ \ i/(9+i-p) 



- - \Ai\y\)pJ 

and 

max{|W(y)| : y e n B3/2 \ B2/3} < Mmax{|?;^(^)| : z S 17^0 B2 \Bi/2} , 
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where M is uniformly bounded because the curvature of d^i is bounded, we obtain that 

3/2 \ ^2/3 



Dv^{y) is uniformly bounded by some constant C on fl B^m \ i?2/3- Because Dv^{y) 



e(<i+i)/{q+i-p)ov{£y), it follows that 

C 

\Dv{x)\ < 



^l/g+l-p |^|(9+l)/(9+l-p) 

By the mean value Theorem, and using the fact that v vanishes on dfl \ {0}, we derive 

, , Cp(x) 

which implies (|2.7 p . □. 

The construction of solutions of the quasilinear equations (|2.1 p with prescribed isolated 
singularity on the boundary of a general bounded domain Vl is settled upon similar 
constructions when the domain is either a half space, or a ball. 

Proposition 2.3 Assume N ~ 1 < q < 2N — 1 and let H = M.^ = {x ~ {xi, xn) '■ xn > 
0)} and k > 0. Then there exists a unique positive solution u = uf^ G C^{H\ {0}) of i2.1 ]) 
in H which vanishes on dH \ {0} and satisfies , 

u{x)^k^{l + o{l)) asx-^Q. (2.10) 

Proof. Since the function x i— *■ fcxjv is A'^-harmonic in H and vanishes on dH \ {0}, it 
is a supcrsolution of (|2.1 p . We write spherical coordinates in under the form 

X = {(r,cr) e [0,c») X S-^-i = (r, singer', cos 0) : a' G 5^"^^ G [0,7r]} , (2.11) 

then ^ 

Du = Ur'l H — Vo-U, 

r 

where i = a;/ Vc denotes the covariant gradient on and equation (|2.1 p takes the 

form 



(Ar-2)/2 



(2.12) 

„.(Ar-2)/2 \ , ^ ' 

-~r-'^div„. [ [ui + I y^wr V<,u + Iwl' m = 0. 



Next ^ 

sm (!) 
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where e is derived from xj \x\ by the rotation with angle 7r/2 in the plane 0,x,N (N being 
the North pole), and Vo-' is the covariant gradient on S^~'^ and (see [3]) 



sin^-^0(^2 + -f + u^l (2.13) 

(^-2)72 



u 



2 iw |2 



\ya'u\ 



{N~2)/2 



sin^ d) \ \ ' sin 



sm 6 \ \ r2 sm 6 / 



If tt depends only on r and 0, (|2.1 p takes the form 

(Af-2)/2 



-r-2sin2-^, 



(^sin^-2 ^ (m2 + r-2u2y'^ +\ur\ = 0. 



(2.14) 



S'iep _Z We look for a local subsolution w under the form 

w(r,cr) = fc(l-r")r"^cos0 r > 0, 0e [0,7r/2]. 
where a > is to be determined. Then 

w,. = — A:r^^(l + (a — l)r") cos (/) and = —kr~^ {1 — r°') sin (f) 
wl + r-^wj := P = k^r-"^ (l + 2(acos2 - l)r" + r2"((a2 _ 3^) cos^ + 1)) 

Wrr — A;r^'^(2 — (a — l){a — 2)r") cos0 and = — fcr~"'^(l — r") cos0 
Pr = -2k^r-^ [2 + (4-Q!)(acos2(/)- l)r" + (2-a)((a^ - 2a) cos^ (/) + l)r^°'] 

= -k^ar""-^ [2 + (a - 2)r"] sin 20, 
PrWr +r-^P4,w^ = 2^^^^^ [2 + (5a- 6 + (2a - a^)cos2 <?!))r" + 0(r2"))] cos0, 

(A^ — l)r^^ti;r + Wrr + [N — 2)r^^ cot </> + r^'^w^^ 

= [4 _ 2N) + (2 - a)(iV + a - 2)r"] cos 

Since 

—div (^Dw\^ ^ -Dw^ + w'^ ~ Lw 

= _p(A'-2)/2 |-(^_ ly-'^w^J^Wrr + [N - 2^-"^ COi (j) + T 

N ~2 
2 

and 



= fc«(l - r")«r"« cos' </> = fc«(l - qr" + 0(r-'"))r-9 cos« 
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a straightforward computation leads to 

Lw = kP-^a [3 - 2A^ + (2 + a){N - 2) cos^ cj) + 0(r")] p(Af-4)/2^a-7 ^ 

+ k'^{l - qr" + 0(r2"))r-9 cos"? (j) 
= kP-^a [3 - 2Af + (2 + a){N ~ 2) cos^ 0] ^-(2^^-!)+" cos + fc«r-9 cos« 
„ gfc9r-9+" cos? + C)(r.-(2JV-i)+2" cos 0) + 0(r-'?+2" cos 0). 



(2.15) 



By assumption q < 2N — 1. If we choose a < min{2Af — 1 — q, 1/{N — 2)}, there exists 
R e (0, 1] such that Lw < on H n Br. 

Step 2 Next we construct a solution ur in Br n H which vanishes on dBji n H and on 
dH \ {0} and satisfies 

l^'J^^M^k. (2.16) 

r-»0 cos 4> 

Let = A:(l — R°')R~^. Inasmuch w — £r is a subsolution, for any e > we can construct 
a nonnegative solution of (|2.1 |) in iJ n (i?i{ \ S^) which vanishes on 7? n OBr and on 
dH n (-B/j \ -Be) and takes the value ke^'^x^ on H O dB^. By comparison 

{w{x) - £r)+ < u,{x) < kxN \x\~'^ . (2.17) 

Furthermore, e i-^ is increasing. Set u — ur = lim^^oUe, then m is a solution of (|2.1 p 
in iJ n Br which vanishes on OBr n H and on dH \ {0} and satisfies the same inequality 
p. 17 p as Ue, but in whole HHBr. This implies that ()2.16 p holds uniformly on [0, n/2-S], 
for any (5 > 0. In order to improve this inequality, we perform a scaling: for r > 0, we set 
u^{x) = ru{rx). Then satisfies 

- div (iT^u'-r^' Du"-^ + r2^-i-9(M'')« = (2.18) 

in iJ n BR/r where there holds 

k{xN \x\'^ (1 - la;]") - £fl.)+ < m'^(x) < kxN \xf^ ■ (2.19) 

Since is uniformly bounded for 1/2 < |a;| < 2, it follows from regularity theory [12] that it 
is also bounded in the C^'"-topology of 2/3 < |a;| < 3/2. Using Ascoli's theorem and the fact 
that u'^{x) converges to kxN pointwise and locally uniformly, it follows that Du^{x) — 
r'^Du{rx) converges uniformly in {x £ H : 2/3 < |x| < 3/2} to —2kxN \x\~* x + k \x\~'^ 
which is the gradient of I . Using the expression of Du in spherical coordinates 

we obtain 

r^Ur'i — rudjG H — — -Va'U — 2fccr7vi + keN uniformly on S!^^^ as r ^ 0, 
sm (j) 

where = (c, ejv). Inasmuch i, e and V(j'U are orthogonal, the component of e^r is sine/), 
thus 

ru^{r, a' ~i —k sin (f> as r ^ 0. (2.20) 
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Since 



u{r,a',(j)) = / u^{r,a\9)d0, (2.21) 

J7r/2 



the previous convergence estimate establishes (|2.16 p . 

Step 3 Construction of the solution in H . Let 77 be the truncation function introduced in 
the proof of Proposition 12.11 and ?/e(r') = rj{{r — e)+). Then the function uj^^^ defined by 
UR^t = in H D Br and zero outside, is a subsolution of (|2.1 p in H which vanishes on 

dH \ {0} and satisfies (|2.16 p . Using the same device as in Step 2, we construct a sequence 
of solutions Us {S > 0) of (|2.1 p in H \ Bg with boundary value k5~'^XN on dBs H H, zero 
on dH \ Bg and satisfies 

UR.e < US < kxN \x\"'^ . 

When 5^0, us decreases and converges to some u which satisfies (|2.1 p and the previous 
inequality. Letting successively e — > and ri{r) wc obtain that u satisfies 

ur{x) <u{x) <kxN\xf^ inH, (2.22) 

where u is the extension of u by zero outside Br . The proof of (|2.10 p is the same as in Step 
2. 

Step 4 Uniqueness. Let u and u be two solutions of (|2.1 p satisfying (|2.10 p and e > 0. Then 
Me = (1 + e)u + e is a super solution which is positive of dH \ {0}. Inasmuch it dominates u 
both in a neighborhood of and in a neighborhood of infinity, it dominates u in H. Letting 
e — + yields to u > u. Similarly u > u. □ 

Proposition 2.4 Assume N - 1 < q < 2N - 1 and let B = Bi{0), a e dB and k > 0. 

Then there exists a unique function u — u^ ^ G C^{B \ {a}) which vanishes on dB \ {a} and 
satisfies h2.1 ]) in B and 

1 - \x\ 

u{x)^k ^-^{l + o{l)) as X a. (2.23) 

\x — a\ 

Proof. With a change of coordinates, we can assume that B has center m ~ (0, 0, —1/2) 
and a is the origin of coordinates. We denote by lo the point (0, 0, —1) and by X^^ the 
inversion with center uj and power 1. By this involutive transformation, the half space 
H = {x ^M.^ : XN > 0} is transformed into the baU S* = {x G : \x\^ + xn < 0}. Thus 
the function x 1— > Pk{x) = —k{\x\^ + xiq)/2 \x'f' is iV-harmonic and positive in B* , vanishes 
on dB* \ {0} and is singular at 0. Let be the solution of (|2.1 p in H satisfying (|2.10 p . 
and Uk — Vk oT^. Then u^ G C(B* \ {0}) satisfies 



-divi^Duu^ ^ Dufc^ + |x — a;| m^ = in i?* 

Wfc = on dB* \ {0}. 



(2.24) 



Furthermore u\^ < Pk and 

Pk{x) = fclZi_^_H! ^ ^ l/2-|x-m| ^^ ^ ^^^^^ ^ ^^^^^^^ ^ ^^^^^ ^2.25) 
2 \x\ 2 \x\ 
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as a; ^ 0. Inasmuch |a; — wj < 1, is a subsolution of (|2.1 p in B* . For e > we construct 
a solution of (|2.1 p in i?* \ -Be(O) with boundary value Pk- By the maximum principle 
Wfc < We < -Pfe in _B* \ Bg{0). Since the sequence {wg} is monotone, we obtain that there 
exists a solution linie^o :— u ^ C^{B* \ {0}) of (|2.1 p in B* which satisfies 

Uk{x) < u{x) < Pk{x) mB*, (2.26) 

and 

y(^) = fci/^— ^^^(l + o(l)). (2.27) 

We change the variables in setting x'j^ — xn + 1/2 and x'^ = Xi [i = 1, ...^N — 1). We 
define u'{x') — u{x) and denote by a the point (0, 0, 1). Clearly u' satisfies (|2.1 p in -B1/2, 
vanishes on dBi/2 \ {a} and 

w'(a;) = fc^/^— M_ (1 + 0(1)) as a/2. (2.28) 
2\x — a/2\ 

By the transformation £ ^ £P/(9+i-P)u'j.(fe), where £ = 1/2, we obtain a solution Uk a of 
(|2.1 p in B which verifies 

Ufc,a(a;) = 2^/fa+i-^)fc ^~ ^"'^ (1 + 0(1)) as a; a. (2.29) 
\x — a\ 

Because k is arbitrary, (|2.23 p follows. Uniqueness of the solution is obtained as in Propo- 
sition with Ug = (1 + e)u. □ 



Proposition 2.5 Assume N-l<q<2N-l and let G = B , a e dB and k > 0. Then 
there exists a unique function u — ^ G C^{G \ {a}) which vanishes on dB \ {a} and 
satisfies h2.1 ]) in G and 

\x\ - 1 

u{x)^k^ 2(1 + 0(1)) as 2: ^ a. (2.30) 

\x — a\ 

Proof. Uniqueness follows from (|2.30 p by the same method as in Proposition l2.3l and Propo- 
sition [T4l Actually, it will be proved in Theorem l2.7l For existence we perform the inversion 
Iq with center and power 1. It transforms the function ^ constructed in the previous 

proposition into a function v G (G \ {a}) which vanishes on dB \ {a} and satisfies ()2.30 p . 
Furthermore v is solution of 

- div (\Dv\^-^ Dv) + l^r'^ \v\''-^ v^Q (2.31) 

in G. Since \x\ > 1, u is a super solution for (|2.1 p in G. With no loss of generality we can 
assume that a = (0, ...0, 1) and let be the solution of (|2.1 p in = {x = {xi, xn '■ 

XN > 1)} satisfying ()2.10 p already constructed in Proposition 12.31 Then — rj{u^ ^) is a 
subsolution in G (where rj^ has been defined in the proof of Proposition 12. ip . By the same 
approximation as in the previous proposition, we construct an increasing sequence {u^} 
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(e > 0) of solutions of (|2.1 p in G \ B^{a) which vanishes on dG \ B^(a), takes the value v 
on G n dB^{a) and verifies < < v m G \ B^{a). Letting e ^ 0, we obtain the existence 
of a solution u* in G which satisfies 

<u* <v in G (2.32) 

where we denote by u^l^ the extension of u^l^ by zero in H^"^ . We conclude that (|2.30 p 
holds in H^. In order to extend this convergence to whole G, we proceed as in the proof 
of Proposition 12.3^ with a minor modification due to the geometry. We put the origin of 
coordinates at a, takes the same spherical coordinates and obtain again that 

r^u*i — rule H Vo-'U* — 2fc(TAri + fee miiformly on S^^^ as r — > 0. 

^ sm (j) 

Therefore (|2.20 p holds for any e [0, 7r/2]. For r > 0, the angle (p ranges from tplr) = 
cos^^(— r/2) to (here is the difference with the half-space case) and \x\ Vm(x) remains 
bounded in this domain, by the regularity theory for quasilinear elliptic equations. Since 

u*ir,a\cl))^ r u;{r,a',9)d0, (2.33) 



we derive, as in the proof of Proposition 12. 3i 

lim u*(r, cr', 0) = fc cos uniformly on [0, 7r/2]. (2.34) 

The proof that (|2.30 P holds is a particular case of Theorem 1 2. 71 □ 

In a general domain we have to extend the solution through the boundary. We denote 
by p(x) the signed distance from x — *■ dVl, that is p{x) = p{x) if x € fl and p{x) = ~p{x) if 
a; € ri"^. Since dVl is G^, there exists /3o > such that if x € verifies — /3o < p{x) < /3q, 
there exists a unique G dfl such that \x — ^x\ — \p{x)\- Furthermore, if is the outward 
unit vector to d^l at ^x, x — — p{x)v(^^ . In particular (^x — pix)i>^^ and ^x + p{x)v^^ have 
the same orthogonal projection £^x onto dft. 

Let Tp„{n) ^ {x <eR^ : -/3o < pi^) < Pa}, then the mapping H : [-Po,l3o] x dO. ^ 
T/3o(f2) defined by Il{p,^) — ( — P^iO is a G^ diffeomorphism. Moreover Dn(0,^)(l,e) = 
e—iy^ for any e belonging to the tangent space T^{dil) to dft at ^. If a; e Tpg{il), we define the 
reflection of x through by tp{x) = £,x+p{x)vxi^ ■ Clearly ip is an involutive diffeomorphism 
from TinT,3„{n) to fl" D T^o(^)- Furthermore for any ^ e dfl, Dip{Cj = Sr^^du) is the 
symmetry with respect to the tangent space T^{dfl) to dfl at ^. If a function v is defined in 
n n T/3o(r2), we define d in nT) T^^ifl) by 

^(x) = <^ (2.35) 
1 -?;oV'(x) if xe rj'=nT/3„(f7). 
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Proposition 2.6 Let v € C'^^°'{n n Tp^iVL) \ {0}) he a solution of IKH) m no Tfj^ifl) 
vanishing on dfl \ {0}. Then v e C^'°'{Tf3g{fl) \ {0}) is solution of a quasilinear equation 



J2^Aj{x, Di) + b{x) |{}|«"^ i = 

OXj 



(2.36) 



in Tp^{VL) \ {0} where the Aj and b are functions defined in T/3g{il) where they verify 

' (i) Aj{x,0) = 



{Hi) ^ 



dAj . , 
or\i 



(2.37) 



<r|?? 



p-2 



. (w) r > 6(a;) > 7 
/or aZZ X € r^(f]) \ {0} /or some /3 G (0,/3o], ?? G M^, ^ G arirf some < 7 < T. 

Proof. The assumptions (|2.37 p implies tliat weak solutions of (|2.36 p are C^'", for some 
a > Q [17] and satisfy the standard a priori estimates. As it is defined the function v is 
clearly in r^o(n) \ {0}. Writing Dv{x) = -D{i} o V'(x)) = -D^I;{x){Dv(tljix))) and 
i = ip{x) — ip~^{x) 



dx. 



But 



7?^(7?«).i?^(i?C) = E E 



Si/'i dv 
dxk dxi 




dxk dxj 



EE 



dij^i dt/jj dv I d( 
dxk dxk dxi j dxj 



We set b{x) = 



A,(x,77) = |DV'||DV'(r?)r'E E 



dxk dxk 



and 



A(a;,ry) = (Ai(a;, ry), A^(x, 77)) = |i?7/^(ry)|^-2 p^)t^^(^) 



(2.38) 



(2.39) 



For any f G 951, the mapping -DV'on(0 the symmetry with respect to the hyperplane 
T^{dfl) tangent to dil at ^, so |D'0(^)| = 1. Inasmuch Dip is continuous, a lengthy but 
standard computation leads to the existence of some f3 G (0, (3o] such that (|2.37 p holds in 
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Tfi{Q)nn . If we define A (resp. b) to be \ri\^ rj (resp 1) on Tfj{Cl)nQ and A (resp. \D^j 
on Ti3{n) n ff, then inequalities (|2.37 p are satisfied in T^(il). □ 

Remark. Notice that, similarly to the p-laplacian, the vector field A is positively homoge- 
neous with exponent p — 1 with respect to rj. Furthermore, if for r > we set A'^{x,ri) ~ 
Aj{rx,ri) , then satisfies the same estimates (|2.37 p as Aj, uniformly in rr-i^(r^^r2), for 
< r < 1. Furthermore 

lim ^^(x,??) = |r?|P"^?7. y-qeM.^, Vj = l,...,iV, 

r— *0 



and this limit is uniform on the bounded subsets of . 



pN 



Theorem 2.7 Let ft be a bounded domain with a boundary and a G dfl. Assume 
N — 1 < q < 2N — 1 and denote by p{x) the distance from x to dfl. Then for any k > 
there exists a unique function u — Uk,a G C{Q \ {a}) which vanishes on dfl \ {a}, is solution 
of i2.1 |) and satisfies 

Uk.a{x)^k ^^^\ (l + o(l)) asx^a. (2.40) 
\x — a\ 

Proof. Uniqueness follows from (j2.40 p by the same technique as in the previous propositions. 
For existence let be a ball of radius R such that B]^ C $7 and a S dB\^, and let loi be 
its center. We denote by the solution of (|2.1 p in B\^, which vanishes on dB]^ \ {a} and 
satisfies 

U'[x)^k?—^^^^{l + o{l)) asx^ a. (2.41) 
|a; — a\ 

If we set Us — ris{U'^), we have already seen that Ug, the extension of Us by zero outside its 
support, is a subsolution of (|2.1 p in il. Because V^^, the A^-harmonic function element of 
C{fl \ {a}) vanishing on dil \ {a}, satisfies 

V^ix) = /^(l + o(l)) as X ^ a, X e B\,, (2.42) 
\x — a\ 

there holds kV^ > Us. If = \ {Be{a)} (e > 0), we construct a solution Ue £ C{fle) of 
(|2.1 p in r^e, which vanishes on dfl \ Be{a) and takes the value kV^ on dBe{a) n ft. By the 
maximum principle e is increasing and Us < < kV^ in Vl^. Letting e — > we obtain 
that Uc converges in the Q^^-topology of \ {a} to a solution u = Uk.,a of (|2.1 p in Vl. It 
follows from the previous inequalities that 

Us{x) <u{x) <kV^{x) Va;en\{a}. (2.43) 



In order to prove the asymptotic behaviour, we proceed as in Proposition 12.41 with the help 
of the reflection principle of Proposition 12.61 We fix the origin of coordinates at a = and 
the normal outward unit vector at a to be — ejy. If ii is the extension of u by reflection 
through dVl, it satisfles 

- ^^ij (a;, Du) + b{x) it = (2.44) 
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in T^{^1) \ {0}. For r > 0, set u^(x) = ru{rx). Then is solution of 

'A 



Y^A'^ix, Du') + r2^-i-96(rx) lil^'"' = (2.45) 



in Tl^-'~\{V) \ {0}, where rj'' := r-^ll. By T, Th 2.4] there exists C > such that 

p{x) 



kV^\x) < Ck'- 



|2 ■ 

a; 



Furthermore, for any x £ T^{n) \ {0}, p{x) := dist (a;,r2) = ^(VXa;)) (we recah that i^{x) 
is the symmetric of x with respect to (957 as it is defined in Proposition 12. 6[) . and c\x\ < 
\ip{x)\ < \x\ for some c > 0, the same relations holds if r^(f7) is replaced by T^'' {n'') 
and p{x) by Pr{x) :— dist (x,r2''). Since O is C^, 

lim - ^ 



1-^0 rpr{x) 

uniformly on bounded subsets of M.^ . Consequently 

l^'^l (,) < Ckr-'B^ = Ck^f^il + 0(1)). 

For < a < 6 fixed and r < ro (for some tq E (0, 1]) the spherical shall Fa.b — {x E : 
a < < is included into T*^'' ^(^1'"). By the classical regularity theory for quasilinear 
equations [17] and Proposition [2?6l there holds 

ll^"nic=(r.,3,3/.)^^^ll"niL==(r.,.,)' (2-46) 

where Cr remains bounded because r < 1. By Ascoli's theorem and (|2.43 p u^{x) converges 
to kxN \x\ in the C^{B^^2 \ ^i/2)-topology. This implies in particular 

lim Du{rx) — —2kxpfx\x\ ^ + k\x\ ^ejy. 

If we take in particular \x\ — 1, we derive 

lim(r{t(r, cr), r^Vt2(r, cr)) — (fc cos i/i, — fc sin (/jgn), (2.47) 

r— >0 

uniformly with respect to cr = (sin(/)CT',cos0) G 5^^^ x [0,7r]. Because 9f2 is there 
exists eo > and a real valued function h defined in Oej, :— B^g n (we recall that 
dH = {x — {x', 0)}) and an open neighborhood Ve^ of such that dV, n Ve^ = {x = {x',xn '■ 
xm — h{x')}, and Dh{0) = (this expresses the fact that dH — To{dil)). If we define ^ by 

'^{x)^{x',XN~h{x')) V.TeV,„. 

then det{D^) = 1 and I?^(0) = /. Up to replacing eo by a smaller quantity, \1/ is a 
diffeomorphism from V^g into a neighborhood V of such that V^g O 9r2) = Q^o- Because 
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dist{'i'{x),dH) = Xn - h{x'), dist = p{x){l +o(l)) as a; ^ 0. Thus, if we set 

X = \E'~-^(y) and u{x) — u*{y), (|2.47 p is equivalent to 

lim i\y\u*{\y\,a),\y\^Vu*{\y\,a)) = (kcos^^-ksin^e^), (2.48) 

uniformly on S^^^, thus 

|2/|u*(|2/|,(7) = fcsin0(l + o(l)) as \y\ ^ (2.49) 

uniformly with respect to cr G S^~^, because u* vanishes on D dH \ {0}. This implies 
(|CTD . ■ □ 

Clearly the mapping k i—t u^^a is increasing. As Uk satisfies the estimates (j2.7 p and 
(|2.30 p . Mfc^a converges in the Cl^^{Q \ {a})-topology, as A: — > oo, to some Moo, a, solution of 
p.l p in n, vanishes on dH. \ {a} and satisfies 

lim \- - '^I'^o^A^) ^ ^. (2.50) 

x^a p[x) 

In order to describe the precise behaviour of Uoo,a, we have to introduce separable solutions 
of (|2.1 p in \ {0}: if we look for solutions u under the form u{r,a) = r^uj{a), then 
P = -(3q = -N/{q +l~ N) and uj satisfies 

(2.51) 

on 5^ ^ where A = (iV — 1)/?^. We shall denote by Sq the set of (always C^'") solutions of 
(|2.51 p . If u is a separable solution of (|2.1 p in iJ which vanishes on dH \ {0}, the function 
a; is a solution of (|2.51 P in S^^-^ which vanishes on dS^^^ = We shall denote by S* 

the set of such functions and by S*^ the subset of positive solutions. We recall some simple 
facts 

Proposition 2.8 (i) For any q > N — 1, Sq contains at least the three constant functions 
anrf±((Ar- l)/3f)i/(9+i-^). 

(li) For any q>2N -\, S* = {0}. 

(Hi) For any q G (A^ — 1, 2 A — 1) . S*j_ contains a unique element. 

Proof. Assertion (i) is evident since A > 0. Assertion (ii), as well as the existence part of 
assertion (iii), can be found in or [52]. Furthermore any w G 5*^ is positive in S^^^ 
and verifies < by Hopf boundary lemma as the outward normal derivative on dS^^^ 
is d /d<j). We can construct a minimal element in 5*^ in the following way: If we denote by 
the unique solution of (|2.1 p in H which satisfies (12.10 p and set Tr{u^){x) — r^^u^ (rx) 
for r > 0, then Tr{u^) is a solution of (|2.1 p in H which satisfies 

r,(uf ) = r(2^-i-«)/(«+i-^)fc^(l + o(l)) as X ^ 0. 
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Thus Tr{u^) = u^^2N-i-q)/{q+i-N)i^- Furthermore, if e + j the maximum principle at 
and at infinity (replacing Uuj by + e and letting e — s- 0) leads to 

w^(r, a) := r-''«cj(a) > uf (r, a) V(r, a) £ (0, oo) x S'^-\ Vfc > 0. 

Letting k ^ oo imphes u^{r,a) > u^{r,a) and Tr{u^) — given that 2N — 1 — q > 0. 
Then the function is invariant with respect to the transformation T^. It is therefore self- 
similar, and consequently under the form u^(r, cr) — r^^''ui{(T). As a result of the previous 
inequality uj is the minimal element of S*_^_. Next we denote S* — max{(5 > : Suj < ui} 
and u^^s* = S*u^. Notice that S* g (0,1] as w > in S^~^ and satisfies Hopf boundary 
lemma on dS!^~^ Clearly u^^S' is a subsolution for (|2.1 p and it is dominated by in H. 
Furthermore S*lu < in S^^^, 5*lo^ < lo_^ on dS^~^ ^ and 

(i) either there exists (Tq G S^^^ such that (5*w((To) = ^ii(o'o), 

(ii) or (5*cj < in S^^^ and there exists CTq G 5*^^^ such that (5*W0((Tq, 7r/2) = tti0(o'o, 7r/2). 

In case (i), and as Du^ never vanishes in iJ, it follows from [BJ Lemma 1.3] (a variant of the 
strong comparison principle) that u^^s* = u- This implies that Ui^^S" is a solution, S* — 1 
and, consequently uj ^ lj. 

In case (ii) we follow the linearization procedure already introduced in 6 . By the mean 
value theorem 

I _f/|-'V-2 I |7V-2 / H \ 

j 

where 

aij = \tiDu^ + (1 - U)Du^j'\^ \UDu^ + (1 - ^O^^c^J:** |^ 

+ (iV - 2) + (1 - tO^^c..*. .0 [Uu^^^ + (1 - , 

with < ti < 1. Next w = — Wcj.i* is positive in H and satisfies 

y 

where c = {{u^Y — ul^ S')li'^oo~'^^,s*) > 0- Notice that (cty (x)) is the Hessian of a strictly 
convex function therefore it is nonnegative and that (a^ )(r, CTq, 7r/2) is positive-definite. 
Therefore it is positive-definite in a neighborhood of (r, (To,7r/2) (independent of r, actu- 
ally). Inasmuch (u^ — u^,S')xn = at (r, ctq, 7r/2), we derive a contradiction with Hopf 
lemma. Therefore case (ii) cannot occur and = w. □ 

Remark. If we look for separable solutions of 

- div {\Du\^'^'^ Duj + lul"^^^ u = 0, (2.52) 
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in R^, where q>p— l>Q,p not necessarily equal to N or to 2, under the form u{r, a) = 
r^u}{a), then /3 = f3p^q — —p/ {q + 1 — p) and uj is a. solution of 

// „N (p-2)/2 \ / „s (p-2)/2 , 

-dzz;, f (/32_,C^2 ^ |y^^|2 j j ^) |'^2 ^^2 ^ |y^^|2 j ^ Q 

(2.53) 

on S^~^ where A{p,q) — Pp^^{q/3p.q — p). If we look for separable solutions in H which 
vanishes on dH \ {0} the solution oj of ()2.53 |) is subject to the boundary condition w — Q 
on dS^~^ = S^^^. A fairly exhaustive theory of existence is developped in [25], 0. The 
existence of non-trivial solution of (|2.53 P is insured as soon A(p,q) > 0, or equivalently 
q < N{p - 1)/{N ~ p) if p < N, and no condition if p > A^. If g > N{p - l)/(iV - p) 
no solution exists, up to the trivial one. This is linked to the removability result proved 
by Vazquez and Veron [18j . The existence of non trivial solutions of the same equation in 
5^^"'^ vanishing on dS!^~^ is much more complicated. However it is proved in [52], [S] that 
there exists a critical exponent qc > p — 1 such that, if q > no non-trivial solution exists 
while a p — 1 < q < qc there exist a unique positive solution in S^~^ vanishing on dS^~^. 
The uniqueness proof in the previous proposition is valid. 

The next result characterizes the solution of ()2.1 p with a strong singularity on the 
boundary. In order to express the result, we assume that the outward normal unit vector to 
dfl at a is — e^r. 

Theorem 2.9 Let il be a bounded domain with a boundary and a € Assume 
< p—1 < q < 2N — 1. Then there exists a unique function u £ C^{fl \ {a}) which vanishes 
on dn \ {a}, is solution of 12.1 ]) in fl and satisfies 

\x - u{x) 

lim \ = oo. (2-54) 

x^a p[x) 

Furthermore 

lim \x-af''u{x)=uj{a), (2.55) 

X — -a 
{x — a) I \x — a\ (7 

locally uniformly on Finally 

Proof. We already know that Uoo,a satisfies (|2.54 p . By translation we fix the origin of 
coordinates at the point a and we assume that — ejv is the outward unit vector to 951 at 

0. If G is any domain in to the boundary of which belongs, we denote by the 
solution of (|2.1 p in G, which vanishes on dG \ {0} and verifies 

uG^/c^^£M(l + o(l)) asa;^0, (2.56) 
\x\ 

where [x) — dist (a;, G). When there is no ambiguity, — u^. By the maximum principle 
G C G" implies < in G. By dilation we can assume that there exist two balls of radius 

1, B C and B' C with respective center b = and 6' = —6 with the property that 
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= OBndB'. It follows from the maximum principle, the fact that (x) = iS{x)) where 
S is the symmetry with respect to the hyperplane dH and Proposition [2131 Proposition [^3] 



i^) u^{x) <Ukix) <ui\x) <u^ b' + 



x-b' 
\x-b'\^ 



S b' 



yx& B 



(m) Ukix) < u^" (x) < ui { S { b' + 



and similarly 



x-b' 



{^) uiix)<uj^{x)<urix)<u^{S{b' + 



x-b' 



yxe B 



(h) uj^ix)<ur{x)<u^l^sl^b' + ^^-^ 
Letting /c — > oo, we obtain 



{i) ugix) < uUx) < ug' (x) < \^b' + ■ 

(ti) Uooix) < u^\x) < ug (^S (^b' + 
as well as 

{^) ugix)<\x\-''^cjii^]<u^'{x)<ug{s{b' + 



x-b' 



yx e B 



Vx e n, 



X 

(u) \x\-''^cj[fj]<u^\x)<ug(s(b' + 



\x-b'\' 



yxeB 



x-b' 
\x-b'f 



Va; e H. 



From (|2.60 p -fi) and the fact that b' = —b, we also derive 



\xr^u{x/\x\)<u^\x)< 



sU±K-b 



\x + b\' 



-A 

LU 



/ S (x + b-\x + b\^ 



But 



S 



\x + b\' 



- b 



\x + b\ 



S[x + b-\x + 6r b 
|a;|(l + o(l)) asx^O 



(remember that |5| = 1). If x = {xi, . . . , xn), \x + b\^ = \x\^ + 1 + 2xn and 



S (^x + b — \x + b\^ b^ — {xi, . . . ,xn + 



(2.57) 



(2.58) 



(2.59) 



(2.60) 



(2.61) 
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Thus ^^IW} becomes 

\x\-^'uj{x/\x\)<u^^ix) < \x\-'^^ \x + bf'u 



X + \X\ Gat 



(2.62) 



x\Jl + \xr + 2XN 



If we assume = °ixN) then (a: + eAr)/(|x| y 1 + \x\^ + 2xn) = x{l + o(l))/ |x| as 
X ^ 0, and 

If we define T by 



uZ\x) = \xr>'^u;{x/\x\)il + o{l)). 



(2.63) 



then (|2.60 reads also as 



T (x) ijj ' 



|r-i(x) 



<u>^{x)<\x\ '''^(A 



(2.64) 



Furthermore 



Then 



r-\x) 



Xl 



and 



\r-\x)\ 



p,,fr-\x) 



\x-hY 



\x-by 



V|T-i(x)| 



\x\~'''' \x-hpuj 



p i X- \x\ ejv 



If we assume again \x\^ = o[xn) then {x — \x^ ei^)/{\x\ y 1 + \x^ — 2xn) = x{l + o(l))/ 
as a; ^ 0, and 

ug{x)^\xr^'^u;{x/\x\)(l + oil)). (2.65) 
Combining (|2.59 ^ -(i). (|2.62 p and (|2.64 p we obtain that 

Moo(x) = |a;r^''tj(a;/|2:|)(l + o(l)) as x ^ (2.66) 

uniformly on any subset of ^l such that \x\^ = o(xjv) near 0. In order to obtain the precise 
behaviour (|2.55 p . we proceed and in the proof of Theorem 12.71 We extend u by reflection 
through d^l near and denote by u the extended function defined in Tf^{il). For re (0, 1] 
we define 

Wr := Tr{u){x) — r^''u{rx). 

Then Wr satisfies 



' dx 



(2.67) 



19 



in Since Wr is uniformly bounded on Ti/2,2 (by Proposition 12 . II applied to u and 

—u) and the definition of the refected function), Dwr{u) is bounded in C"(r2/3^3/2)- By 
Ascoli's theorem Wr converges in the C"'^(r2/3^3/2)-topology to x 1-^ |a;|~'''' a)(a;/ |a;|), where 
uj is defined from lo by reflection through the equator dS^~^. In ordre to get rid of the 
boundary, we use again the diffeomorphism ^ which sends B^^ onto itself and verifies 
*(Beo n dn) = n oh. We set X = *^H2/) and u{x) = u*{y). Then 

lim {\yf-u*{\y\,a), |y|^'+^ Vu*(|y| , a)) = (^(0), -c.^ gn), (2.68) 

uniformly on S^^^, thus 

u*i\y\,a)^\yf-Lomi + o{l)) as \y\ ^ 

uniformly with respect to cr G S^~^, because u* vanishes on H dH \ {0} 
stronger result than (12.55 p follows, namely 

u(a;) = |a;p'^«w(a;/|x|)(l + ol)) as a; 0. (2.70) 

Mutatis mutandis, this estimate implies uniqueness of a solution with a strong singularity 
as in Theorem 12.71 □ 



(2.69) 
Actually, a 



3 The removability result 

In this section SI is a domain of and a £ dil. The next result extends the removability 
result of Gmira-Veron dealing with semilinear equations. 

Theorem 3.1 Let g be a continuous function defined on R which satisfies 

liminf g(r)/r'?'= > and limsupg(r)/|r|'^'^ < 0, (3.1) 

where qc '■= 2N — 1 and let u € C^{fl \ {a}) be a solution of 

~ div (^\Duf^^ Du^ + g{u) = inn (3.2) 

which coincides with some (f) £ C^{drt) on dfl \ {a}. Then u extends to fl as a continuous 
function. 

Proof. Without any loss of generality, we can assume that SI is bounded, a = and — e^v is 
the outward normal vector to dil at 0. We denote by the solution of (12.2 p in SI which 
vanishes on 9S1 \ {0} and satisfles 

n"(^) = r$(l + °(l)) asx-.0. 

Let M be the supremum of \4>\ on i9Sl and M = max{Af, (B/A)^^'}. By assumption there 
exists A > and B > 0, depending only on g, such that 

-div(\Duf^'^Du)+Au''-'<B m{xen:u{x)>0}. (3.3) 
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liv = u-M, then v <0 ondn\ {0} and 

- div (\Dvf^^ DvJ + Av''" < m{xefl: v{x) > 0}, 



(3.4) 



Using the same functions i]^ as in the proof of Proposition [27T] we deduce that rj^iv) satisfies 
the same inequahty as v, but on whole fl. By Proposition l2.2l with q = Qc and the expression 



of Vp it follows that 



(3.5) 



where the constant c depends on A and N. Furthermore, there exists a function u* G 
Ci(n\ {0}) such that < w+ < u*{x) < cV^ in n, and 



div ID 



,N-2 



Du* ] + Au* ''■^ = inn. 



(3.6) 



As in the proof of Theorem 12.91 we extend u* through the boundary into u and scale it by 
setting Tr{u) := Wr{x) = ru{rx) for r > 0. Inasmuch all the previous a priori estimates 
apply (compactness), it follows that there exists a subsequence {r„} converging to and a 
function w £ C'i(M^ \ {0}) such that Wr^ — !■ w in the Cj^^^-topology of \ {0}, w is a 
solution of 



(1^ 



-div ( \Dw 



Dw] + Awl" ^0 in \ {0} 

w>0 in H ^ {x eR'^ : XN > 0} 
w = ondH\{0}. 



(3.7) 



At end, p. 5 p transforms into 



< w(x) < Vx £ H. 

\x\ 



(3.8) 



For e > we denote by the solution of 

-div {\DW^\^~'^ DW^ + AW^" =0 inH\ 5,(0) 

W, = ce'^XN on HC] dB,{Q) 
W^=(} on OH \B^{0). 



(3.9) 



By the maximum principle < w{x) < Wc(x) < cxm \x\ for any e > 0, and by uniqueness, 
Tr{W(){x) = rWt{rx) = W^/r{x). Furthermore e ^ is increasing. Letting e — > we 
conclude that decreases to some Wq, which is a solution of 



—div ^1 



AW, 



la _ 



Q inH 



Wo>Q in H 
Wa^Q ondH\ {0}, 



(3.10) 
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by the standard regularity results, and satisfies < w < Wq- Finally, Wq inherits the 
following scaling invariance property Tr{Wo){x) = Wo{x) for any r > 0. Therefore Wq is a 
separable solution which endows the following form 

Woix)^Wo{r,a)^r-'u,ia), 

where lu is nonnegative on S^~^ and satisfies 



// n\(N-2)/2 \ , „s ( 



m b^ 



= on dS^-^. 

(3.11) 

By Proposition [JUl = 0. Thus VFo = w = 0, which implies Wr(x) Q a& r ^ 'id and 
equivalently ru{rx) in the C/^^-topology of \ {0}. Consequently Du(x) = o[\x\ ") 
as a; — > and finally u*{x) — o{V^{x)) as x — s- 0. The maximum principle and the positivity 
of u* yields to u* = and finally u < M in In the same way u > —M. Because the 
modulus of continuity of u is uniformly bounded near 0, by the classical regularity theory 
of degenerate elliptic equations (see |12| for example) , u extends as a continuous function in 
whole SI. □ 



4 The classification theorem 

The next result extends some of Gmira-Veron's classification theorem [7l Sect. 4, 5 ] obtained 
in the study of problem (jl.3 p . In the above mentioned article, the main idea was to reduce 
the equation to a infinite dimensional quasi-autonomous evolution system in x S^~^ and 
to use Lyapounov-energy function. Such an approach cannot be adapted in the quasilinear 
case. Our method is based upon scaling and uniqueness arguments. 

Theorem 4.1 Assume N — 1 < q < 2N ~ 1, is a bounded domain with a boundary, 
a G do. and — ejv is the outward normal unit vector to dil at a. Let u e C^(f2 \ {a}) be 
a positive function satisfying 112.1 \) in f2 and vanishing on dfl \ {a}. Then the following 
alternative holds. 

(i) Either there exists k > such that 

u{x)^k ^^^^ , (l + o(l)) asx^a. (4.12) 

\x — aj" 

Furthermore u — Uk,a, the unique solution of h2.1 \) defined in Theorem \2. 7[ 

(ii) Or 

u{x) = \x~a\~'^'' Lu{a){l + o{l)) as x ^ a. (4.13) 

where lo is the unique positive solution of h2.51 \] on S^^^ which vanishes on dS^^^, in 
which case u — Uoo.a- 
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Proof. We assume that a = with z/q = — bn and define 

u(x) ,. u(x) 
J, ^ = lim sup sup — fT— , 

x^O V^'{x) r^O \x\=7- V^\x) 



fc = hmsup ^^^^ ^ = hmsup sup ttq-t^- (4.14) 



Suppose fc = 0. It follows from the maximum principle that for any e > there exists a 
sequence r„ ^ such that < u{x) < €V^{x) in \ {i?r„(0)}. This fact implies the nullity 
of u. Therefore we assume that /c ^ 0. Assume first that k is finite. Then, for any e > 0, 
there exists a sequence of points a;„ converging to such that 



and 



lim -^^^ = k (4.15) 



u 



(x) 



™P TT^ <k + e. (4.16) 
Since Uk satisfies (|2.40 p with a = 0, the two previous relations can be replaced by 



W hm„^oo — 7 — r = 1 

Uk[Xn) 

u( x) 

Hi) s^PN<.„ ^ < 1 



(4.17) 



We denote r„ = |a;„|, ^„ = a;„/r„ and define u„ = rnuirnx) and Ukn ^ r„Uk{rnx). By 
the previous arguments combining a priori estimate and regularity theory, there exist a 
subsequence {rnj} and two nonnegative functions v and u', A^-harmonic in H and vanishing 
on dH \ {0}, such that {unj,Uknj) converges to {v,v') in the C/Qj,-topology of H — M^. 
Clearly equality (|2.40 p implies that r„^. Uk [ruj x) converges to kV^ (x) (which is defined by 
kV^ [x) := kxN / \x\ ^) in the same topology. Since v' is uniquely determined by its blow-up 
at 0, this implies v' = fcVjj^ in H. Furthermore there exists ^ € S^'"^ such that ^ ^. 

If e e ^+'\ «(e) = «'(?), while, if C e 91-/91/(0 = 9,„«(0 = ai;V9i/(e). In 

both situation, the tangency conditions of the graphs of v and v' and the strong maximum 
principle implies that v — v' — kVf^ . By estimate (14.17 l) -(i) and the convergence properties, 
it follows 

lim — ^ "^'^ = 1, uniformly on 1^1 — 1. 

Consequently, for any S > 0, there holds, 

(1 - S)uk{x) < u{x) < (1 + S)uk(x) \/x en\Br„, 

for n large enough, which leads to Uk = u. At end we consider the case fc = oo. Writting 
(|2.1 p under the form 

- di^; (li:*?/!^"^ Dm) + dix)u^-^ = (4.18) 

where d{x) = \u\'^^^ ^ {x) < Cjxp^, by (|2.4 p . We use the boundary Harnack principle. 
By [21 Th 2.2] there exists a constant c ~ c{N, q,il) > such that 

1 u(y) u(x) u(y) / , , „s 

— < ——- < c-^ (4.19) 

c p(y) p{x) p{y) 
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for any a; and y in such that = \y\ be small enough. Since there exists a sequence 
Xn ^ such that lim„^tx) u{xn)/V^{x„-j oo, this implies that 

lim Jinf^^^ : = |a;„|l = c5o. (4.20) 



V"{x) 

Thus u satisfies ((^371) : Theorem [gH and ([^351 imply that holds. □ 

The assumption of positivity on u can be weakened if a better a priori estimate is already 
known. The next result extends [B] Th 1.2] into the framework of boundary singularities. 

Theorem 4.2 Assume N — 1 < q < 2N — 1, il is a bounded domain with a boundary, 
a G dfl and — e^v is the outward normal unit vector to dil at a. Let u ^ C^{il\ {a}) be a 
solution of V2.1 \) in f2 vanishing on dVL \ {a} such that u/V^ is bounded in fl. Then there 
exists fc g R such that u — u^^a- 

Proof. The outline of the proof are very similar to the finite case of the previous theorem. 
We still assume a = and define k by (|4.14 p . If fc = the maximum principle implies u < 
and we return to Theorem 14. II in the case u<0. Iffc^^O, fc>0 for example, (|4.15 p and 
(|4.16 p apply. By the previous scaling method we derive that converges to some function 
V in the C/^^-topology oi H — which is TV-harmonic in H and vanishes on dH \ {0}. 
Because rn^.Uk{rnkX) converges to kVf^ , the tangency condition of v and kVf^ at some ^ 
implies that v = kVf^ . Thus u{x) > for |a;| — r„^. for Uk large enough. This implies that 
M > in and we are back to Theorem 14. II □ 

Remark. In the semilinear case of problem (|1.3 p , it is proved in [7J that any signed solution 
u which satisfies lima;_,a |a; — u{x) = has constant sign. The exponent N characterize 
the minimal changing sign harmonic function vanishing on dQ, with an isolated singularity 
at a. Changing sign singular A^-harmonic functions are constructed in [3 . In particular 
there exist singular iV-harmonic functions w under the form 



w{r, a) = r ^^uj{a) 



where 



_ 7iV - 1 + V]V2+T2]VTT2 
" 6(iV- 1) 



and Lu is defined on ^ = {x e ^ : xn > 0}, vanishes on the equator dS^ ^, is 
positive on S^~^ n {x : xn~i > 0} and negative on S'^~^ n {x : xn-i < 0}. A natural 
question is therefore wether any signed solution u of ()2.1 p in fi which vanishes on dft \ {a} 
and satisfies liuirc^a \x — a\^^ u{x) = has constant sign, and can be henceforth classified 
through Theorem 14. II 

Final remark. If one replaces the A^-harmonic operator by the p-harmonic operator (p > 1) 
and tries to extend the results of sections 2, 3, 4, several difficulties will appear. Even if the 
existence of separable singular solutions is known, the precise value of the exponent /3 > 
such that (r, a) i— > r~'^(p{a) is p-harmonic and positive in H and vanishes on dH \ {0} is 
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unknown but for the specific cases iV = 2 or p = iV or p = 2. Notice that in that case the 
function (f> satisfies the so-called spherical p-harmonic spectral equation 



-div„ 



I / n\ (p-2)/2 \ / „\ (p-2)/2 



= on as'+ . 

(4.2f) 

where A = - 1) + p - A^). lip = 2 then fi = N -\, while if A?^ = 2, /? is the positive 

root of the equation 

3/3^ + 2^^—4/3 -1 = 0. (4.22) 
p- 1 

Furthermore, up to now and due to the lack of conformal invariance, it has not been possible 
to construct the equivalent of the in a general smooth bounded domain that are 
positive p-harmonic functions in il, vanishing on dSl \ {a} and satisfying 

lim \x- af u{x) = <^{a). (4.23) 

X — > a 
(x — a) I \x — a\ — > (J 

However, if 17 = = M;'^ the removability and the classification results of Sections 3 and 4 
are still valid. The proofs of these theorems are developed in [T]. 
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